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Abstract 

Context. Stellar-mass black holes offer what is perhaps the best scenario to test theories of gravity in the strong-field 
regime. In particular, f(R) theories, which have been widely discus in a cosmological context, can be constrained 
through realistic astrophysical models of phenomena around black holes. 

Aims. We aim at building radiative models of thin accretion disks for both Schwarzschild and Kerr black holes in f(R) 
gravity. 

Methods. We study particle motion in /(7?)-Schwarzschild and Kerr space-times. 

Results. We present the spectral energy distribution of the accretion disk around constant Ricci scalar f(R) black holes, 
and constrain specific f(R) prescriptions using features of these systems. 

Conclusions. A precise determination of both the spin and accretion rate onto black holes along with X-ray observations 
of their thermal spectrum might allow to identify deviations of gravity from General Relativity. We use recent data on 
the high-mass X-ray binary Cygnus X-l to restrict the values of the parameters of a class of f(R) models. 
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1. Introduction 

General Relativity (GR) is consistent, in some cases with 
good precision, with observational results (see for instance 
Will 2006). However, an energy momentum tensor repre- 
senting exotic matter (loosely called 'dark energy', e.g. Li, 
Li, and Wang 2011) must be introduced in the right hand 
side of Einstein's equations to fit the currently available 
data when these are interpreted in the framework of the 
standard cosmological model (based on GR). Dark energy 
can be modelled by a cosmological constant, or by a scalar 
field with an equation of state given by p — ludeP, where 
w DE < -1/3 (Biswas et al. 2010a, 2010b). None of these 
descriptions is free of problems, since the energy density 
associated with the cosmological constant that is inferred 
from astronomical observations is approximately 120 orders 
of magnitude lower than the value predicted by field the- 
ory (e.g. Weinberg 1989, Capozzicllo and Faraoni 2010), 
whereas the scalar field has features that are at odds with 
the scalar fields of particle physics (Sotiriou and Faraoni 
2010). 

A different approach to explain the cosmological data is 
to modify the field equations of the gravitational field, in 
such a way that the ensuing theory differs from Einstein's 
in the low-curvature regime. Since there is no a priori fun- 
damental reason to restrict the gravitational Lagrangian 
to a linear function of the Ricci scalar R (see for instance 
Magnano et al. 1987), more general theories can be formu- 
lated using nonlinear functions of this scalar. The so-called 



f(R) theories (e.g. Capozziello and Faraoni 2010) were first 
used to mimic the inferred accelerated expansion of the uni- 
verse by Capoziello (2002). Currently there is in the litera- 
ture a handful of /(-R)-models in agreement with available 
data (De Felice and Tsujikawa 2010). 

Although the present revival of f(R) theories is mainly 
due to their use in the description of phenomena that take 
place for low values of the Riemann curvature, these theo- 
ries have also been applied to gravity in the opposite regime. 
As there is no direct evidence of the behaviour of the grav- 
itational field for very large values of the curvature, the 
early universe and compact objects offer the possibility 
to find deviations from GR. Among the studies in modi- 
fied gravity in the strong regime, we can mention the suc- 
cessful inflationary model based on the R + aR? theory 
(Starobinsky 1980), and the related studies of reheating 
(Motohashi and Nishizawa 2012) and particle production 
(Arbuzova, Dolgov, and Rcverbcri 2012) in the early uni- 
verse. Also of importance is the treatment of neutron stars 
(Cooney, DeDeo and Psaltis, 2010) and black hole solutions. 

Different aspects of black hole physics in f(R) theo- 
ries have been discussed in the literature by Psaltis et al. 
(2008), Hendi & Momcni (2011), Myung (2011), Myung et 
al. (2011), Moon, Myung, and Son (2011a, b), Hendi et al. 
(2012), and Habib Mazharimousavi et al. (2012). Static and 
spherically symmetric black hole solutions were obtained 
via perturbation theory by de la Cruz-Dombriz, Dobado, 
and Maroto (2009), whereas black holes with these symme- 
tries have been studied by means of a near-horizon analysis 
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by Perez Bergliaffa & De Oliveira (2011). Finally, f(R)- 
Kerr-Newman black holes solutions with constant Ricci 
scalar have been recently studied by Ccmbranos and col- 
laborators (Cembranos et al. 2011). 

From an astrophysical point of view, thin accretion for 
the Schwarszchild space-time in /(i?)-gravity has been dis- 
cussed by Pun et al. (2008), without including the expected 
spectra of concrete astrophysical black holes or any com- 
parison with observational data. 

In the present work we investigate the existence of stable 
circular orbits in Schwarzschild and Kerr f(R) space-times 
with constant Ricci scalar and analyze the main features 
of accretion disks around these black holes. In particular, 
we present temperature and spectral energy distributions 
for Page-Thorne disks, and compare the results with those 
obtained using the standard Shakura-Sunyacv model. 

The paper is organized as follows. In Section 2 we pro- 
vide a brief review of f(R) theories of gravity. Circular or- 
bits in both Schwarzschild and Kerr f(R) space-times with 
constant Ricci scalar are studied in Section 3. Section 4 
is devoted to the calculation of the properties of accretion 
disks in these space-times. Consequences for some specific 
prescriptions for the function / are discussed in Section 5. 
We close with some considerations on the potential of as- 
tronomical observations to test the strong-field regime of 
gravity. 

2. f(R) gravity 

In f(R) gravity, the Lagrangian of the Hilbert-Einstein ac- 
tion, given by: 



is generalized to: 
%] 



16ttG 



16ttG 



(i? + /(i?))V^d 4 x, 



(1) 



(2) 



where g is the determinant of the metric tensor, and f(R) 
is an arbitrary function of the Ricci scalar. In the metric 
formalism the field equations are obtained varying Eq. (2) 
with respect to the metric: 



R^{l + f\R))-\g^{R + f(R)) 



(3) 



16ttG 



r„„ = o, 



□ = V^, f(R) = 
df(R)/dR, and the energy momentum tensor is defined by: 



where i? M „ is the Ricci Tensor 



(4) 



Here, C m stands for the matter Lagrangian. 

Equations (3) are a system of fourth-order nonlinear 
equations for the metric tensor g^ v . An important difference 
between them and the Einstein field equations is that in 
f(R) theories the Ricci scalar R and the trace T of the 
energy momentum tensor arc differentially linked, as can 
be seen by taking the trace of Eq. (3), which yields: 



R(l + f'(R)) - 2(R + f(R)) - 3D/'(i?) - 



16ttG 



T = 0. (5) 



Hence, depending on the form of the function /, there may 
be solutions with traceless energy-momentum tensor and 
nonzero Ricci scalar. This is precisely the case of black hole 
space-times in the absence of a matter source. 

Notice that in the case of constant Ricci scalar Ro with- 
out matter sources, Eqs. (3) can be re-written as: 



where: 



and, by Eq.(5): 



A 



Ra 



/(go) 
2f(R ) 



(6) 



(7) 



/W-i 

Hence, in this case, any f(R) theory is formally, but not 
physically, equivalent to GR with a cosmological constant 
given by Eq.(6) . 

We shall investigate the existence of stable circular or- 
bits in Schwarzschild and Kerr f(R) space-times with con- 
stant Ricci scalar in the next section. 

3. Circular orbits around a black hole in f(R) 
3.1. f(R)-Schwarzschild space-time 

The Schwarzschild space-time metric in f(R) theories with 
constant Ricci scalar Rq takes the form (Cembranos et al. 
2011): 



ds 1 = - 



, 2GA/\ 
cf- 



c 2 R 



12 



2 

r 



(8) 



dr 2 



[(1- 



Ro r 2] 
12 ' J 



where Rq, given by Eq. (7), can take, in principle, any real 
value. Since we are looking for space-time metrics that may 
represent astrophysical black holes, we shall select those 
values of Ro that lead to acceptable solutions. 

The radius ro of the horizon follows from the condition 
<7oo( r o) = 0- From Eq. (8), the values of r corresponding to 
the event horizon satisfy: 



c 2 R r 3 - 12c 2 r + 24GM = 0. 

In terms of the following adimensional quantities: 
r 

x = — , 
Ro = i?or|, 

where r g = GM/c 2 , this equation takes the form: 
R x 3 -12x+24 = 0. 



(9) 

(10) 
(11) 

(12) 



We show in Figure 1 the Ricci scalar as a function of the 
radial coordinate of the event horizon. We see that for Ro <G 
(0,4/9) there is an inner black hole event horizon and an 
outer cosmological horizon, whereas for Ro < there is only 



1 Accretion through thick disks onto Schwarzschild and Kerr 
black holes with a repulsive cosmological constant was studied 
by Rezzolla, Zanotti and Font (2003), and Slany and Stuchh'k 
(2005), respectively. These studies do not analyze spectra or 
compare results with observational data. 
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Figure 1. Plot of the Ricci scalar as a function of the ra- 
dial coordinate of the event horizon in /(i?)-Schwarzschild 
space-time. 



one black hole event horizon. The event and cosmological 
horizons collapse for Ro = 4/9 and for larger values of the 
Ricci scalar naked singularities occur and hence, no black 
holes are possible 2 . We shall, then, restrict ourselves to the 
study of trajectories for Ro £ (-co, 4/9]. 

3.1.1. Equations of motion and effective potential in 
/(i?)-Schwarzschild space-time 

The geodesic equations in the metric given in Eq. (8) can be 
obtained by means of the Euler-Lagrange equations using 
the Lagrangian: 



L 



2GM\ c 2 R 



r j 



12 

■2 



(13) 



2GM' 



+ r 2 (6 2 + sin6 2 ^ 



where x^ = dx^ /da, and a is an affinc parameter along the 
geodesic x fl (a). The resulting geodesic equations for t and 
4> are (setting = it/ 2): 



2GM\ R 



12 



c 2 r J 



r 2 $ 



(14) 
(15) 



where k and h are constants. An equation for r that is 
simpler than the one obtained from the Lagrangian fol- 
lows from the modulus of the 4-momentum p, given by 



9pv 



x^x 



e 2 , where e 2 



for massive particles, and 



2 The results obtained in our analysis of /(i?)-Schwarzschild 
space-time with constant Ricci scalar are consistent with those 
given by Stuchh'k and Hledfk (1999) in Schwarzschild-de Sitter 
and Schwarzschild-anti de Sitter space-times. 



e 2 = for photons. It takes the form: 



2 2GM\ c 2 R 2 



r J 

^2 



12 



(16) 



[(1 



2GM^ 



Eo r 2] 
12 1 J 



2 12 2 

r = e , 



with i and 4> given by Eqs. (14) and (15), respectively. The 
set of Eqs. (14), (15), and (16) completely determine the 
motion of a particle in the /(i?)-Schwarzschild space-time. 

3.1.2. Trajectories of massive particles 

Equations (14), (15), and (16) can be used to obtain the 
so-called energy equation (e.g. Hobson et al. 2006): 



r 2 + 



2GM 



2GM\ 
c 2 r J 



c 2 R 

i 

12 



Ro ■ 

r 

12 

c 2 (k 2 



(17) 



!)■ 



The constant k is defined as k — E/itiqc 2 , where E repre- 
sents the total energy of the particle in its orbit, and rooc 2 
its rest mass energy. The constant h stands for the angular 
momentum of the particle per unit mass. From Eq. (17) we 
can identify the effective potential per unit mass as: 



Feff (r) = — (l - — 
2r 2 \ c 2 r 



_v 

12 



(18) 



2GM c 2 R 



12 



■0 2 

r 



The extrema of the effective potential are obtained by look- 
ing for the roots of the derivative of the latter equation with 
respect to the radial coordinate. In terms of x, Ro, and the 
adimensional angular momentum per unit mass of the par- 
ticle h = /i(cr-g) -1 , this reads: 



eff 



dx 



v3 



3h^ 

x 4 



1 

+ -- 



Rpx 
12 



= 0. 



The derivative of Eq. (19) with respect to x gives: 

d 2 V eS _ c 2 (-4R x 4 + 15R x 3 + 12x - 72) 
^k 2- ~ ~i? x 4 (-| + 1) ' 

where we have replaced h by (Harko et al. 2009): 



(19) 



(20) 



(21) 



where x c corresponds to the radius of a circular orbit. We 
have performed the numerical calculation of the extrema of 
the effective potential for different values of Ro- 

For Ro > 0, as shown by Stuchlik ct al. (1999) and 
Rezzolla et al. (2003) in Schwarzschild-de Sitter space-time, 
stable circular orbits exist for values of the specific angular 
momentum that satisfy: 





_ Ro 


v ^ 


12 




_ 3_ 





< h < h 



(22) 



where hi sco and h osco stand for the local minimum and local 
maximum of the specific angular momentum at the inner 
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and outer marginally stable radii. From Eq. (20) we see 
that the existence and location of the circular orbits depend 
on the Ricci scalar Ro. By equating Eq. (19) to zero and 
isolating Ro, we obtain the Ricci scalar as a function of the 
radial coordinate of the circular orbits: 



12 (6 - x c ) 
(15 - 4x c ) x c 3 ■ 



(23) 



In Figure 2 we show the plot of the latter equation. Wc 
see that there is an upper limit for the Ricci scalar, Ro = 
2.85 x 10 -3 , for which circular orbits are possible. In Figure 
3 we plot the effective potential as a function of the ra- 
dial coordinate. The dots indicate the location of the in- 
nermost stable circular orbits. The corresponding values of 
the event and cosmological horizons, radii of the innermost 
and outermost stable circular orbits, for six different values 
of R G (0, 2.85 x 1CT 3 ) arc shown in Table 1. Wc sec that 
for increasing values of the Ricci scalar the event horizon 
becomes larger than in Schwarzschild space-time in GR as 
well as the location of the innermost stable circular orbit. 



0.004 r 



0.003 - 



p| 0-002 



0.001 - 



0.000 




Figure 2. Plot of the function given by Eq. (23). The 
absolute maximum corresponds to x = 15/2 and Ro = 
2.85 x 10~ 3 . 



The extrema of the effective potential for Ro < are all 
located outside the event horizon, as shown in Figure 4. The 
value of the radial coordinate for the event horizon is less 
than 2 (i.e. smaller than for Schwarzschild black holes in 
GR). The location of the innermost stable circular orbit is 
closer to the horizon than that of the Schwarzschild solution 
in Einstein's gravity. The limit Rq — > — oo in Eq. (20) yields: 



x 3 (-4x+15) = => x = V x = 3.75. 



(24) 



The plot of the effective potential corresponding to the val- 
ues of the parameters of Table 2 is shown in Figure 5. 



3.2. Kerr space-time in f(R) theories 

The axisymmetric, stationary and constant Ricci scalar 
geometry that describes a black hole with mass, electric 
charge and angular momentum was found by Carter (1973), 
and was used to study f(R) black holes by Ccmbranos et 
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Figure 3. Effective potential for different values of Ro > 
and h in /(i?)-Schwarzschild space-time. The dots indicate 
the location of the innermost stable circular orbit. 
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Figure 4. Plot of Ricci scalar as a function of the radial co- 
ordinate of the event horizon (line) and of the Ricci scalar 
as a function of the radial coordinate of the innermost sta- 
ble circular orbits (dashed line) for Ro G [—1.5,0.45] in 
/(-R)-Schwarzschild space-time. 



al. (2011). The form of the metric is the following: 



ds 2 = f dr 2 



A r Ag 



c dt 



(25) 



Ag sin 2 6 
P 2 

A r f cdt 



(r + a 

2 
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Table 1. Location of the event and cosmological horizon, and of the innermost and outermost stable circular orbits for 
Ro > in /(i?)-Schwarzschild space-time. Here x = r/r g . 



Ro 


Radius 
event horizon 


Radius 
cosmological horizon 


Radius innermost 
stable circular orbit 


Radius outermost 
stable circular orbit 


10" w 


Xeh = 2 


x ch = 3.46 x 10 b 


Xisco — 6 


Xosco = 14421.70 


10" 6 


x c h = 2 


x ch = 3463.10 


Xi S co = 6.00016 


Xosco = 143.45 


10" 4 


Xeh = 2 


Xch = 345.40 


xisco = 6.02 


Xosco = 30.16 


ten 3 


x oh = 2.00067 


Xch = 108.53 


xisco = 6.19 


Xosco — 13. IT 


2 x 10 -3 


Xe h = 2.00134 


Xch = 76.44 


xisco = 6.51 


^osco — 9.80 


2.84 x KT 3 


x eh = 2.0019 


Xch = 63.98 


xisco = 7.40 


Xosco — T.60 



Table 2. Location of the event horizon, and of the inner- 
most stable circular orbit for Ro < in /(i?)-Schwarzschild 
space-time. Here x = r/r g . 



Ro 


Radius 
event horizon 


Radius innermost 
stable circular orbit 


-i<r a 


Xch = 1-999 


xisco = 5.86 


-10" 2 


Xch = 1.993 


Xisco = 5.26 


-ht 1 


x ch = 1.939 


Xisco = 4.35 


-1.5 


Xch = 1.541 


Xisco — 3.83 



14r 

12 
10 



> 6 





Ro=-l(T 3 h=3.4940 




Ro=-l(T 2 h=3.7050 




Ro=-10 -1 h=4.8612 




Ro=-1.5 h=l 1.9075 
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Figure 5. Effective potential for different values of Ro < 
and h in /(i?)-Schwarzschild space-time. The dots indicate 
the location of the innermost stable circular orbit. 



where: 



1 -Ro 

1 — — — r 



12 



2GMr 



p 2 = r 2 + a 2 cos9 2 , 
A e = 1 + -y|a 2 cos 2 i 



5 = 1 



Ro 
12 



a 2 . 



(26) 
(27) 
(28) 

(29) 



Here M and a denote the mass and angular momentum per 
unit mass of the black hole, respectively, and Ro is given 
by Eq.(7). 

Because of the constancy of p 1 and p^ along the trajec- 
tories, and of the reflection-symmetry of the metric through 
the equatorial plane, the orbit of any particle with initial 
condition p e = will remain in the plane tt/2, where the 
metric has the form: 

„2 _2 



ds 2 = 



2 )dt 2 



r 



dr 2 



(30) 



2ac 

AA2 



- A r ) dtdcj) 



Here, 



S = 1 



Ro 
~V2~' 



2GMr 



R 



2 

12 a • 



(31) 
(32) 



If Rn 



0, Eq. (30) represents the Kerr space-time metric 
in GR as expected. 

The equation that yields the position of the event hori- 
zon is obtained by setting l/g IT = 0: 

2GMr „ 



r' + a" 



1 - -7-r 



R 2 
r 



In terms of x, Ro, and a = 
form: 

(x 2 +a 2 ) (l 



12 

Rpx 2 
12 



0. 



this equation takes the 



2x = 0. 



(34) 



In Figure 6 we plotted the Ricci scalar as a function of the 
radial coordinate of the event horizon for Ro € [—0.3, 1] 
and a = 0.99 (i.e. a nearly maximally rotating black hole, 
such as Cygnus XI, Gou et al. 2011). If R € (0,0.6], there 
are 3 event horizons: the inner and outer horizons of the 
black hole and a cosmological horizon; for Ro > 0.6 there 
is a cosmological horizon that becomes smaller for larger 
values of Ro. If Ro € (—0.13,0) there are 2 event horizons. 
For Ro < —0.13 naked singularities occur. In the following 
we shall analyse the existence of stable circular orbits for 
R G (-0.13,0.6]. 

3.2.1. Equations of motion and effective potential in 
/(i?)-Kerr space-time 

In order to obtain an expression for the effective potential, 
we make use of the invariant length of the 4-momentum p: 

gTPvfr, = e 2 , (35) 
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Figure 6. Plot of the Ricci scalar as a function of the 
radial coordinate of the event horizon for Ro € [—0.3, 1] 
and a = 0.99 in /(i?)-Kcrr space-time. 

where e 2 = c 2 for massive particles and e 2 — for photons. 
Since we are only interested in trajectories on the equatorial 
plane, we set pe = 0, and Eq. (35) gives: 



g u (p t ) 2 + 2^ + 9* 
where: 
Pt = kc 2 , 

Pr Qrrf • 



2 , rr 2 2 
6 + .9 Pr = £ , 



(36) 

(37) 
(38) 
(39) 

Substituing Eqs. (37), (38), and (39) into (36), the equation 
for f 2 takes the form: 



9 



- g tt ( fcc 2) 2 + 2g ^hkc 2 -g++h s 



(40) 



The contravariant components of the space-time metric 
given by Eq. (30) are: 



— 2 



9 



A r c 2 r 2 
A r 



A r a 2 



■9"" = ~-2 



cA r r 2 

—2 



if 



A r r 2 



a (r 2 + a 2 - A r ) , 



(A r - a 2 ) 



(41) 
(42) 
(43) 
(44) 



Hence, the energy equation for a massive particle is given 
by: 



\f + V eS (r, a, R , k,h) = - (k 2 - l) 



(45) 



where the effective potential is: 



V eS (r,a, R ,k, h) 



—2 



2r 2 



+ T( fc2 - 1 )-^ ( 46 ) 



and: 



r = [(r 2 + a 2 ) ck -ah} 2 - A r {ack - hf 



(47) 



3.2.2. Equatorial circular orbits of massive particles 

If a massive particle is moving in a circular orbit of radius 
rj SCO , the value of the effective potential at any point of the 
orbit satisfies the equation: 



V e g (r isco , a, Ro, k, h) = — (k 2 - l) 



so. 



dVcs 
dr 



(r,a,R ,k,h) 



For a stable orbit, the equation: 
d 2 V cSi 



dr 2 



(r,a,R ,k,h) 



0. 



>0, 



(48) 
(49) 

(50) 



must also be satisfied. As shown in Section 3.2, there are 
black holes if R € (-0.13,0.6]. 

Numerical calculations of the radius of the innermost 
stable circular orbit for several values of Ro < 0, with a = 
0.99, show that the stable circular orbits lay outside the 
event horizon (see Table 3). Notice that as the value of the 
scalar decreases, the radius of the innermost stable circular 
orbit becomes smaller. The radius of the innermost stable 
circular orbit in Kerr space-time in GR (rj sco = 1.4545 r g , 
in the prograde case) is always larger than in /(i?)-Kcrr. 
In Figure 7, the effective potential that correspond to the 
values of Table 3 are shown. 

We follow Stuchlfk and Slany (2004) to study the exis- 
tence of stable circular orbits for Ro G (0, 0.6]. The specific 
angular momentum of a massive particle in a co-rotating 
circular orbit yields (Stuchlfk and Slany, 2004): 



2a + ax c (x c 2 + a 2 ) § - x c (x c 2 + a 2 ) ( £ 



1/2 



l _ 3 _ a^Ro + 2a I ' 

X c 12 



1/2- 



1/2 



(51) 

We see from the latter equation that circular orbits must 
satisfy the following two conditions: 



X c < 



1/3 



(52) 



which is the same for /(i?)-Schwarzschild space-time with 
positive Ricci scalar, and: 



1 

x c 



12 



2a 



1 



Ro 

12 



1/2 



> 0. 



(53) 



The minimum and maximum of Eq. (51) give the values 
of the specific angular momentum that correspond to the 
innermost and outermost stable circular orbit respectively, 
once the angular momentum a of the black hole and Ro 
are fixed. We show in Table 4 the values of such radii for 
different values of Ro, and a = 0.99, and in Figure 8 the cor- 
responding plot of the effective potential. As expected, the 
radius of inner circular orbit is larger than in Kerr space- 
time in GR. We also found, by equating to zero the deriva- 
tive of Eq. (51), that stable circular orbits only exist if 
R € (0,1.45 x 10- 1 ). 

The analysis of the circular orbits presented in this sec- 
tion will be applied next to the construction of accretion 
disks around black holes. 
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Table 4. Location of the event and cosmological horizon, and of the innermost and outermost stable circular orbits for 
Ro > and a = 0.99 in f(R)-Keri space-time. Here x = r/r g . 



Ro 




Radii 
event horizons 


Radius 
cosmological horizon 


Radius innermost 
stable circular orbit 


Radius outermost 
stable circular orbit 


lO" 


Xchl 


= 0.86 Xch2 = 1.14 


x c h = 3463 


Xi S co = 1.4545 


Xosco = 143.59 


10" 4 


Xehl 


= 0.86 x ch2 = 1.14 


x ch = 345.40 


x isco = 1.4547 


Xosco — 30.53 


6.67 x 1(T 4 


Xchl 


= 0.86 x oh2 = 1.14 


= 113.12 


x isC o = 1.4560 


Xosco = 16.002 


1(T 3 


Xehl 


= 0.86 x ch2 = 1.14 


x ch = 108.53 


xisco = 1.4567 


Xosco — 13.93 


10~ 2 


Xehl 


= 0.85 x oh2 = 1.15 


x c h = 33.59 


xisco = 1.4765 


Xosco — 6.25 


10" 1 


Xchl 


= 0.83 x oh2 = 1.22 


x ch = 9.80 


Xisco — 1.92 


^osco — 3.22 



Table 3. Radii of event horizons and circular orbits for a 
/(i?)-Kerr black hole of angular momentum a = 0.99, for 
some values of Rq < 0. Here x= r/r g . 



Ro 



Radii event horizons 



Radius innermost 
stable circular orbits 



-10~ a 


Xehl 


= 0.86, X e h2 


= 1.14 


Xisco 


= 1.452 


-1.2 x 10" 3 


Xehl 


= 0.86, x c h2 


= 1.14 


Xisco 


= 1.451 


-10" 2 


Xehl 


= 0.86, x ch 2 


= 1.13 


Xisco 


= 1.43 


-to- 1 


Xehl 


= 0.91, X c h2 


= 1.03 


Xisco 


= 1.20 


-1.25 x 10 _1 


Xchl 


= 0.96, x c h2 


= 0.98 


Xisco 


= 1.04 



-0.2 



-0.3 



-0.4 



-0.5 



-0.6 




Ro=-10 -3 h=1.5677 
Ro = -1.2 10" 3 11=1.5674 
Ro^-lO" 2 h=1.5611 
Ro=-10 -1 h=1.4237 
Ro=-1.25 10"' 11=1.2208 



0.6 0. 



1.0 1.2 1.4 1.6 
r/r g 



Figure 7. Effective potential as a function of the radial co- 
ordinate (Ro < 0, a = 0.99), in /(i?)-Kerr space-time. The 
dots indicate the location of the innermost stable circular 
orbit. 



4. Accretion disks in strong gravity 

4.1. Standard disk model in general relativity 

The first realistic model of accretion disks around black 
holes was formulated by Shakura and Sunyaev (1973). They 
considered that the matter rotating in circular Keplerian 
orbits around the compact object loses angular momentum 
because of the friction between adjacent layers and spirals 
inwards. In the process gravitational energy is released, the 



SB 

> 



-0.220 
-0.225 
-0.230 
-0.235 
-0.240 
-0.245 
-0.250 




Ro=10-" h= 1.5683 
Ro=l(T 4 h= 1.5684 
Ro=6.67 l(T 4 h= 1.5689 
Ro= 1CT 3 h= 1.5692 
Ro=10 -2 h= 1.5760 
Ro^lO -1 h= 1.6076 



2 3 4 5 

r/r g 

Figure 8. Effective potential as a function of the radial co- 
ordinate (Ro > 0, a = 0.99), in /(i?)-Kerr space-time. The 
dots indicate the location of the innermost stable circular 
orbit. 



kinetic energy of the plasma increases and the disk heats 
up, emitting thermal energy. 

Novikov, Thorne, and Page (Novikov & Thorne 1973; 
Page & Thorne 1974) made a relativistic analysis of the 
structure of an accretion disk around a black hole. They as- 
sumed the background space-time geometry to be station- 
ary, axially-symmetric, asymptotically flat, and reflection- 
symmetric with respect to the equatorial plane. They also 
postulated that the central plane of the disk coincides with 
the equatorial plane of the black hole. This assumption en- 
tails that the metric coefficients gu, gt<f>, 9rr, gee, and g^ 
depend only on the radial coordinate r. 

The disk is supposed to be in a quasi-steady state 
(Novikov & Thorne 1973), so any relevant quantity (for 
example the density or the temperature of the gas) is av- 
eraged over 2tt, a proper radial distance of order 2H 3 , and 
the time interval At; that the gas takes to move inward 
a distance 2H . During At, the changes in the space-time 
geometry are negligible. 



3 Here H represents a particular height above the central plane 
of the disk (\z\ < H « r). 
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Particles move in the equatorial plane in nearly geodesic 
orbits; consequently, the gravitational forces exerted by the 
black hole completely dominate over the radial accelera- 
tions due to pressure gradients. 

The expression of the energy flux for a relativistic ac- 
cretion disk takes the form (Novikov & Thorne 1973; Page 
& Thorne 1974): 




where Mg stands for the mass accretion rate, fi for the 
angular velocity and E and L represent the specific energy 
and angular momentum, respectively. The lower limit of the 
integral ri SCO corresponds to the location of the innermost 
stable circular orbit. 

The angular velocity f2, the specific energy E, and the 
specific angular momentum L of the particles moving in 
circular orbits are given by (Harko et al. 2009): 



i ■ r 




0.0 I I V i I i I i I 

5 10 15 20 




Figure 9. Energy flux as a function of the radial coordi- 
nate of an accretion disk around a Schwarzschild black hole 
in Shakura-Sunyacv (SS) and Page-Thorne (PT) models, 
respectively. 



E 
L 



d<t> ~9t<p,r + \J (~gt<j>,rf - gtt,r£ 
dt 9<l><l>,r 

9tt + gt<t>& 

\/-gtt - 2^fi - g^fl 2 
9tij> + 9<j><j>Q 

\l~9tt - igt^ - g<f><t>Q 2 



(55) 
(56) 
(57) 



Eqs. (56), and (57) can be derived by writing the effective 
potential V e s (r) in terms of the metric coefficients and solv- 
ing for E and L the equations V e g(r) = and V e g^(r) = 0. 
The formula for the angular velocity f2 = d(j>/dt is obtained 
by substituing E, and L into the geodesic equations dt/ dr 
and dfy/dr (Harko et al. 2009). 

In the next subsections we calculate the energy flux, 
temperature and luminosity of an accretion disk around 
a Schwarzschild and a Kerr black hole in GR and f(R) 
gravity with constant Ricci scalar, adopting the following 
values for the relevant parameters: M — 14.8M©, M = 
0.472 x 10 19 gs _1 , and a = 0.99, which are the best estimates 
available for the well-known galactic black hole Cygnus-Xl 
(Orosz et al. 2011, Gou et al. 2011). 



4.1.1. Relativistic accretion disk around Schwarzschild and 
Kerr black holes 

In order to obtain an expression of the energy flux and 
temperature of the disk, for the Schwarzschild black hole, 
we calculate the angular velocity fi, the specific energy E 
and angular momentum L of the particles in the disk, using 
the metric (Schwarzschild 1916): 



ds 2 



-c 2 { 1 



2GM \ 
c 2 r J 



dr 



i-^VV 

c l r J 



(58) 



r 2 sin e 2 dd> 2 . 



From Eqs. (55), (56), and (57) we get: 



n = 

E = 
L = 



GM 



(1 



2GM S 



3GAV 



1 



3GM ' 



(59) 
(60) 

(61) 



By replacing these equations in Eq. (54), we obtain: 
3M c 6 1 



Q(x) 



SttxV 2 (GM) 2 



1-5 



(62) 



Vx+ \/3tanh 1 



where Xj SCO = 6r g is the location of the innermost stable 
circular orbit in Schwarzschild space-time. 
By means of Stcfan-Bolzmann's law, 



T(r) 



Qir) 

0"SB 



1/4 



(63) 



(where z stands for the correction due to the gravitational 
redshift), the temperature of the disk can be obtained as a 
function of the radial coordinate r. 

In Figures 9 and 10 we plot the energy flux and the 
temperature as a function of the radial coordinate for a 
Schwarzschild black hole in both Keplerian and relativistic 
accretion disk models. We see that GR effects introduce a 
decrease of the peak of the energy flux by a factor w 2, and 
that the temperature distribution is also diminished. 
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5.0x10 s 




Figure 10. Temperature as a function of the radial coordi- Figure 11. Energy flux as function of the radial coordinate 
nate of an accretion disk around a Schwarzschild black hole for an accretion disk around a Kerr black hole of angular 
in SS and PT models, respectively. momentum a = 0.99 in the PT model. 



The Kerr space-time metric (Kerr 1963) in Boycr- 
Lindquist coordinates for 9 = n/2 is: 



ds 2 = -4 (A r - a 2 ) dt 2 + ^dr 2 



2ac 



(64) 



• 2 i 2 

r + a 



+ 



where: 



[r 2 4- 



A r ) dtdcj) 
2 ^~ - A r a 2 



(r 2 + a 2 ) 



2GMr 



(65) 



The expression of the energy flux now becomes: 

M n, r 



Q(r) 



where: 



47r v / ^g 



j (S- QLj L, r dr, 



E - ML) L, r 
M fi jr 



c [(pi +p 2 ) Pg] 
2 



P4 



3M c 2 



P5 



47T A 



E - VLL 



8tt ( Pl + P2 y 



and the coefficients pi are given by: 



Pi = x 

P2 

P3 = (l 



5/2 



<3/2 



,3/2 



a 

x^T 2 



x 3/2 

8a 3a 2 



P4 = X 3 1 - 



2a 



X X 

2 



3/2 



-1 



= -i fx 3 + a - 2ax 3 / 2 ") f 1 - - + 4^r 



x 3/2 



x 3/2 



(66) 
(67) 

(68) 

(69) 
(70) 
(71) 
(72) 
(73) 




Figure 12. Temperature as function of the radial coordi- 
nate of an accretion disk around a Kerr black hole of angu- 
lar momentum a = 0.99 in the PT model. 



Here x = r/r g is an adimcnsional radial coordinate, and 
a = a/r g is the angular momentum of the black hole in 
adimensional units. 

In Figures 11 and 12 we show the plots of the energy flux 
and temperature of an accretion disk around a Kerr black 
hole of angular momentum a = 0.99, whose innermost sta- 
ble circular orbit is at r; sco = 1.4545 r g . In Figure 13 we 
show the luminosity of relativistic accretion disks around 
both Schwarzschild and Kerr black holes. For comparison, 
we also present the Schwarzschild/Shakura-Sunyaev lumi- 
nosity. 

The values of the maximum temperature, luminosity, 
and the energy of the peak of the emission for all these 
models are shown in Table 5. As expected, the highest lu- 
minosity corresponds to a prograding accretion disk around 
a black hole. Since the last stable circular orbit is located 
closer to the black hole than in Schwarzschild space-time, 
thermal radiation is emitted at higher energies. 
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Table 5. Values of the energy of the peak of the emission, the maximum temperature, and luminosity of an accretion 
disk around Schwarzschild and Kerr black holes (a = 0.99) in the SS and PT models. 





Schwarzschild (SS) 


Schwarzschild (PT) 


Kerr (PT) 


^max 


1007.59 eV 


746.94 eV 


1654.9 eV 


Tin ax 


0.406 keV 


0.277 keV 


0.539 keV 




2.42 x 10 37 ergs" 1 


1.22 x 10 37 ergs" 1 


2.26 x 10 37 ergs" 1 



CD 



■ ■ ■ • Schwarzschild (SS) 
— Schwarzschild (PT) 
Kerr(PT) 




2.00 2.25 2.50 2.75 3.00 3.25 3.50 3.75 4.00 

Log 10 (E[eV]) 

Figure 13. Luminosity as function of the energy for a 
relativistic accretion disk around a Schwarzschild and a 
Kerr black hole (a = 0.99). We also plot the luminosity 
as a function of the energy of an accretion disk around a 
Schwarzschild black hole that corresponds to the SS model. 



4.2. f(R)-gravity 

4.2.1. /(i?)-Schwarzschild black holes 

The energy flux of an accretion disk around a f(R)- 
Schwarzschild black hole with metric given by Eq. (8) takes 
the form: 



9A/ c 6 



(i-r 1 



^(GM)\ 5 /36_ 3c 2 R 



E-Q.L) L tI dx 



where 



(74) 



E-flL) L r = - 



y/3 (-12x+72 + 4R x 4 - 15R x 3 



12x 3 



if-Rr 



,!/2 



(75) 

We adopt for the radius of the outer edge of the disk (Dove 
et al. 1997): 

r out = llr isco . (76) 

According to the latter equation, if we take for the inner- 
most stable circular orbit n sco = 6.3 r g , the outer egde of 
the disk yields approximately 70 r g . If larger values of r out 
are considered, there are no major differences in the tem- 
perature and luminosity distributions. 

We first compute the temperature and luminosity spec- 
tra distributions for Ro < 0, and adopt the values given in 
Table 2. In Figures 14 and 15 we show the plots of the tem- 
perature as a function of the radial coordinate, and of the 
luminosity as a function of the energy, respectively. Notice 



that the corrections due to the gravitational redshift have 
been taken into account. 

The maximum temperature as well as the luminosity 
increase for smaller values of Ro (Figures 14 and 15). In 
the four cases presented, the accretion disk is hotter than 
in GR, e.g. in a /(i?)-Schwarzschild disk for Ro = —1.5, 
the maximum temperature and luminosity are a factor 1.9 
and 3.7, respectively, higher. The energy corresponding to 
the peak of the emission shifts to higher values, reaching 
1359.20 eV for the set of adopted parameters. 



5x10 
4x10" 





1 1 I 1 1 1 I 


i | i | i | i 

Ro = -1.5 • 






Ro = - 10" 1 - 




' N \ 
' N \ 
i *v 


Ro = - 10" 2 

Ro = - 10 3 - 
Ro = 




' V 
1 ' ' ■ *v 






1 • if ~"-"-.r;; 






1 ■ i/ 
• / 

: ; 

; j 

i-r . i i i 





10 15 20 25 30 35 40 

r/r n 



Figure 14. Plot of the temperature as a function of the 
radial coordinate for some values of Ro < 0, for a f(R)- 
Schwarzschild black hole. 



We showed in Section 3.1.2 that for Ro > 0, stable cir- 
cular orbits are possible within a minimum and maximum 
radius. We see in Table 1 that only for Ro = 10" 12 and 
Ro = 10" 6 accretion disks arc possible, if we take for the ra- 
dius of the outer edge of the disk r ou t = 70r g . The values of 
the location of the innermost stable circular orbit, location 
in the radial coordinate of the maximum temperature, max- 
imum temperature and luminosity, and the energy of the 
peak of the emission are displayed in Table 7. We conclude 
that for Ro € (0, 10~ 6 ] the temperature and energy distri- 
bution have no significant differences with Schwarzschild's 
distributions in GR. 

4.2.2. /(i?)-Kerr black holes 

We calculate next the energy flux of an accretion disk 
around a /(i?)-Kerr black hole: 



Q(x) 



Mo 



n. 



E - flL) L x dx, (77) 
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Table 6. Values of the location of the last stable circular orbit, location in the radial coordinate of the maximum 
temperature, maximum temperature and luminosity, and the energy of the peak of the emission for an accretion disk 
around a /(i?)-Schwarzschild black hole with Rq < 0. 



/(fl)-Schwarzschild R = R = -10" 3 R = -10"^ R = -KT 1 R = -L5 



fisco/^g 


6 


5.85 


5.26 


4.35 


3.82 


f'TmaxA'g 


10.82 


11.04 


10.05 


7.58 


6.45 


T m ax 


0.277 keV 


0.298 keV 


0.371 keV 


0.468 keV 


0.526 keV 


^max 


746.94 eV 


825.31 eV 


1007.59 eV 


1230.13 eV 


1359.20 eV 


i(-E'max) 


1.22 x 10 37 ergs _1 


1.97 x 10 37 ergs" 1 


2.91 x 10 37 ergs" 1 


3.93 x 10 37 ergs" 1 


4.5 x 10 37 ergs" 1 



Table 7. Values of the location of the innermost stable circular orbit, location in the radial coordinate of the maximum 
temperature, maximum temperature and luminosity, and the energy of the peak of the emission for an accretion disk 
around a /(i?)-Schwarzschild black hole with Rq > 0. 



/(-B)-Schwarzschild R = R = 10~ b 



^"isco / fg 


6 


6.00016 


fTmax / r g 


10.81 


10.82 


-Mnax 


0.277 keV 


0.277 keV 


-E'max 


746.94 eV 


746.94 eV 




1.22 x 10 37 ergs" 1 


1.22 x 10 37 ergs" 1 



and 



37.75 
37.50 
37.25 
37.00 
36.75 
36.50 



Ro = -1.5 


1 1 1 1 1 1 1 1 


; Ro = -10" 1 




_ Ro = - 10" 2 X 


\ \ 


■ Ro = -10"V^'' 


• \ \ 
\ \ 


" Ro = /> 


\ \ 




\ '. \ \ 


/ '/ ■ • ' 


'\ '• * \ 


// .■ 




/// 


\ v ■. \ \ 


/'■' 


\\ ■•■ A ; 


fi' 




Pi 


w ■ \\ 


/' ■' 

f. 1 


\ \ '■ \\ 


ti 1 i.i.i. 


\ \ '. 1 \ 

\\ • \\ 

1 . 1 , \ 1 . . 1 



2.00 2.25 2.50 2.75 3.00 3.25 3.50 3.75 4.00 

Log 10 (E [eV]) 

Figure 15. Plot of the luminosity as a function of the en- 
ergy for some values of Ro < 0, for a /(i?)-Schwarzschild 
black hole. 



where: 



n. x = -36\/3— o?7, 



(78) 



L* = - 



_ x 1 / 2 {l2x 3 + a [12a - aR x 3 - 4V36X 3 -3x 6 R J } 
(12x 3 + a 2 R x 3 - 12a 2 ) 2 a/^FV 3 + 12 

(12 + a 2 Ro)xV7Tni' v ' 

4x [I2x 3 + a 2 (-12 + R x 3 )] (Z 5 + l 6 + h + h + h) 



(12 + a 2 R ) Vl2x 3 - R x«(Zi + hi 



.3/2 



(e - QLj 



12c^/hTT 4 



(12 + a 2 R ) [12x 3 + a 2 (-12 + x 3 R )] 



(80) 
, (81) 



h 

h 
h 
h 

h 
h 
h 
k 
h 

ha 

ho 
hi 

hla 



-72a 3 x - 216ax 3 



12xV36x 3 - 3x 6 R , 



= a 2 V36x 3 



-432a^x 



2„2 



3x 6 R [a 2 xR + 24 + x(l2 + x 2 R )] , 
-x 6 (12 + a 2 R ) 2 + 48a 3 ^36x 3 - 3x 6 R , 



144ax 2 V3fo 



12x 3 [36x 2 + a 4 R 



3x 6 R 

a 2 (36-3x 2 R )] , 



[l08ax 2 (-24 - 12x + 5R x 3 )] ^I2x 3 - R x 6 , 

36a 3 (-48 - 36x + 7R x 3 ) ^/l2x 3 - R x 6 , 

\/3a 4 x (-12 + R x 3 ) [-108 + R x 2 (-3 + R x 3 )] 

36 v^x 5 {72 + x [-12 + R x 2 (-15 + 4x)] } , 

3V3a 2 x 2 {864 + xl 9a } , 
[2160 + xZ 9ol ], 

(432 + R 2 x 4 (15 + 8x) - 12R x (21 + 26x)) , 
144 (-3 + x) x 5 + a 4 R x 3 (-12 + R x 3 ) 
48a 3 \/36x 3 - 3R x 6 , 



144ax 2 ^36x 3 - 3R x 6 



12a 2 x 2 / llQ , 



-36 



-36 + R x 2 (3 + 2x)] 



If we adopt for the radius of the inner edge of the disk 
1.4545r g , according to Eq. (76), r out sa 16r g . From Eq. (77) 
we numerically calculate the temperature and luminosity 
for the values shown in Table 3, taking into account the 
corrections coming from the gravitational rcdshift. The re- 
sults are displayed in Figures 16, 17, and Table 8. We see 
that the temperature of the disk increases for smaller val- 
ues of Ro- The ratio of the maximum temperature between 
the GR and f(R) cases, with R = -1.25 x 10"\ is 1.20. 
The peak of the emission rises a factor of 2, and the corre- 
sponding energy is shifted towards higher energies. 

Since the radius of the outer edge of the disk is 16r g , 
we see from Table 4 that accretion disks are only possible 
for R = 10" 6 and R = 10" 4 , until up R = 6.67 x 10" 4 . 
For such values of Rq , we show in Table 9 the values of the 
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Figure 16. Plot of the temperature as a function of the 
radial coordinate for some values of Ro < of a /(-R)-Kerr 
black hole of angular momentum a = 0.99, corrected by 
gravitational redshift. 



Figure 18. Plot of the temperature as a function of the 
radial coordinate for some values of Ro > 0, for a f(R)- 
Kerr black hole of angular momentum a = 0.99. 



37.400 




32 34 36 Log 10 (E[eV]) 
Log 10 (E[eV]) 



Figure 19. Plot of the luminosity as a function of the en- 
Figure 17. Plot of the luminosity as a function of the en- e rgy for some values of R > 0, for a /(i?)-Kerr black hole 
ergy for some values of R < 0, for a /(_R)-Kerr black hole f angular momentum a = 0.99. 
of angular momentum a = 0.99. 



location of the last stable circular orbit, maximum tempera- 
ture, luminosity, and the energy of the peak of the emission, 
and in Figures 18 and 19, the temperature and luminosity 
distributions respectively. As in /(i?)-Schwarzschild black 
holes with positive Ricci scalar, these differences are minor. 

We proceed now to examine some specific forms of the 
function /, and the constraints imposed on them by the 
previous analysis. 

5. Limits on specific prescriptions for f(R) 

As discussed in Sections 3.1.2 and 3.2.2, the existence of 
Page-Thorne disks around f(R) black holes imposes the 
following limits on Rq: 

— /(i?)-Schwarzschild space-time: 

R e (-oo; 10- 6 ], (82) 



— /(i?)-Kerr space-time: 

R G [-1.2 x 10" 3 ; 6.67 x 10" 4 ]. (83) 

As we shall see in Sect. 6, contemporary observa- 
tions of Cygnus X-l rule out accretion disks around f(R)- 
Schwarzschild black holes, since the maximum temperature 
obtained in such models is lower than the inferred through 
observations (Gou et al. 2011). Hence, we will only con- 
sider the values of Ro given by expression (83). We shall 
show in this section how these values lead to limits on the 
parameters of two examples of f(R) theories via Eq. (7). 
We shall also impose the following viability conditions, to 
be satisfied by any f(R) (Cembranos et al. 2011): 

- 1 < f'(Ro) < 0, (84) 
f"(Ro) > 0. (85) 
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Table 8. Location of the last stable circular orbit and maximum temperature, maximum temperature, luminosity, and 
energy of the peak of the emission for an accretion disk around a /(i?)-Kerr black hole with Rq < and a = 0.99. 



J 1 / n\ T.' ...... 


Ko = U 


D 1 A — 3 

Ko = — 1U 


Ko = —1.2 X 1U 


f isco / fg 


1.4545 


1.4523 


1.4518 


rTmax/r g 


3.79 


3.79 


3.79 


T 

-*- max 


u.uoy K.t. v 




fl "idl 48 keV 
u.tj^i^o ive v 


-E'max 


1659.4 eV 


1659.4 eV 


1659.4 eV 




2.26 x lO^ergs" 1 


2.38 x lO^ergs" 1 


2.41 x 10 37 ergs _1 


f(R)-Kerr 


Ro = -HT* 


Ro = -HT 1 


Ro = -1.25 x 10" 1 


f isco / Tg 


1.4325 


1.2017 


1.0419 


T^Trnax / V* g 


3.85 


3.85 


3.78 


Tux ax 


0.553 keV 


0.663 keV 


0.652 keV 


^raax 


1833.52 eV 


2025.9 eV 


2025.9 eV 


-£'(-£' max) 


2.94 x 10 37 ergs _1 


4.23 x 10 37 ergs _1 


4.60 x 10 37 ergs _1 



Table 9. Values of the location of the last stable circular orbit, location in the radial coordinate of the maximum 
temperature, maximum temperature and luminosity, and the energy of the peak of the emission for an accretion disk 
around a /(i?)-Kerr black hole with Rq > and a = 0.99. 



f(R)-Ken 


Ro = 


Ro = 10~ 4 


Ro = 6.67 x 10" 4 


Tisco J fg 


1.4545 


1.4547 


1.4559 


rTmax/^g 


3.79 


3.79 


3.79 


T m ax 


0.53942 keV 


0.53927 keV 


0.53843 keV 


^raax 


1659.4 eV 


1659.4 eV 


1659.4 eV 


^(^max) 


2.26 x 10 37 ergs" 1 


2.25 x 10 37 ergs" 1 


2.09 x 10 37 erg s" 1 



5.1. }{R)=aRt } 

The parameters a, P and the Ricci scalar are related by Eq. 
(7) as follows: 



Ro 



1 



a 08 -2) 



(86) 



Introducing the adimcnsional parameter a' = RE 1 a with 



Ro 



this equation can be written as 



Ro = 



1 



a' (/3 - 2) 



(87) 



Notice the condition /3 > to ensure the GR limit for 
small values of the Ricci scalar R. Let us consider first the 
case of a positive Ricci scalar, which leads to: 



Case I 



or 



Case II 



a' > A 13 > 2, 
a' < A /? < 2, 
a' < A fi > 2, 



or 



a' > A P < 2, 
and 1/(P — 1) an even number, that is: 



(89) 
(90) 
(91) 

(92) 



with n G Z. By isolating a' from Eq. (87), we obtain the 
function: 



We show in Figure 20 the plot of a' as a function of p (with 
P > 0) for fixed values of the Ricci scalar. We see that for 
P G (0, 2), a' G (-00, 0). For P = 2, Eq. (93) is not defined, 
while for P > 2, a' takes large positive values. Since a needs 
to be small in order to recover GR for small values of the 
Ricci scalar, the case a' > 0, P > 2 is discarded. Hence, we 
obtain the following restrictions on the parameters: 



a'e(-oo;0) A /3g(0;2) A R G (0; 6.67 x 10" 4 ]. 



For negative values of the Ricci scalar, from Eq. (93) we 
require that: 



1-P= 2m, => P odd = 1 - 2m, 



(94) 



1 - P = 2m +1, =S> /3 eV en = -2m, 



(95) 

where m G Zq , so that P > 0. If m = 0, P = 1 and a' = 
a = —1. These values of the parameters lead to f(R) = —R, 
which does not reduce to GR. For /? = 2, Eq. (93) is not 
defined. If p > 3 and is an odd number, a' takes positive 
large values, while if p > 4 and is an even number a' is large 
and negative. Since a' has to be small to recover GR for 
small values of the Ricci scalar, we conclude that negative 
values of Ro are not allowed in this theory. 

We now restrict the values of a and P according to Eqs. 
(84) and (85). The first and second derivative for the given 
/ function are: 



f'(R) = apRP-\ 
f"(R) - aP(p~l)R^ 2 . 



(96) 
(97) 



The restrictions over a and P that satisfy Eq. (85) are: 

a > A P>1, (98) 



a<0 A Pe(0,l). 



(99) 
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R (0;6.67 x 10 ; 



e G (-00; 0) A a' G (Oje' 1 R ), 



or 



0.0 0.2 0.4 0.6 



1.2 1.4 



ee(0;oo) A a' G (e ^ojoo), 
- R G [-1.2 x 1(T 3 ;0) 

e G (-00; 0) A a G (-e" 1 |R | ;0), 

or 

ee(0;oo) A a' G (-00; -e" 1 |R |), 

The first and second derivative of the function / take the 
form: 



P 

Figure 20. Plot of a' as a function of (3 for different values 
of R . 



f'(R) 



1+ln- 

a 



f"(R) = 



(103) 
(104) 



The condition given by Eq. (98) is discarded because it does Th e condition f"{R ) > is satisfied if e > A Rq > 0, or 
not satisfy Eq. (93). The viability condition expressed by e < A R < 0. Equation (84) in adimensional form is: 
Eq. (84) takes the form: 



1 < afiR^- 1 < 0. 



(100) 



1< e 1 + In ■ 



We can constrain the values of a using the latter inequality 
as follows: 

<0 < 1. 

By multiplying by i?o /3_1 01 the latter restrictions yields: 



< 



This equation, together with (104), yields: 



(105) 



> aPRo* 3 - 2 > aRo 13 ' 1 . 
In order to satisfy Eq. (100): 

-1 



(101) 



aRo 



p-i 



> -1 =S> a > 



Ro 



0-1' 



- For e > and R > 0, a' G (e R ; R exp{(± + l)}). 

- For e < and R < 0, 

a' e (-e|R |;-|R |exp{l-^}). 

Summarizing all the constraints, we have that: 

R G (0;6.67 x 10~ 4 ], e > 0, a' G (e _1 Rq; 00), (106) 



If /3 = 0, a > - Ro, and for = 1, a > -1. The set of G f e R o; Rocxp | ^ + 1 
values for a is a G {—Ro, 0). 

We conclude that the values of a and /3 that are permit- , 
ted by our model as well as by the two viability conditions 



are: 



ae(-i? o ;0) A /3G(0;1) A Rq G (0; 6.67 x 10 -4 ]. 



R G [-1.2 x 10~ 3 ;0), e<0, a' G (-e _1 |R |; 0),(107) 



a' G 



-e|R |;-|Ro|exp^l- — 



5.2. f(R) = R eln^ 



The first group of constraints is fulfilled for e > 0, while in 
the second group for only e G (—1/2, 0). We conclude that 
In this case, the parameters e and a, and the Ricci scalar for the f(R) under scrutiny, the values 
are related by Eq. (7) in the simple form: 



a = Rq exp I 1 



(102) 



R G (0;6.67 x 10" 4 ], e > 0, a'e(e R ;oo), (108) 



and, 



Dividing by R g , we obtain a' = Ro exp ( 1) . For R > 0, 

the function a (e) is always positive, while it is neg ativefor Ro G [— 1.2x 10 -3 ; 0), eG (-1/2,0), a' G (-e" 1 |R |; 0) 



all e if Ro < 0. The constrains over e and a that follow from 
Eq. (7) are: 



are allowed. 



(109) 
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6. Discussion 

The results presented in Section 4 can be compared with 
current observational data to derive some constraints on a 
given f(R) theory. In order to illustrate this assertion we 
shall consider Cygnus X-l, which is the most intensively 
studied black hole binary system in the Galaxy. A series of 
recent high-quality papers (Rcid ct al. 2011, Orosz et al. 
2011, Gou et al. 2011) have provided an unprecedented set 
of accurate measurements of the distance, the black hole 
mass, spin parameter a, and the orbital inclination of this 
source. This opens the possibility to constrain modified the- 
ories of gravity with rather local precision observations of 
astrophysical objects in the Galaxy. 

Cygnus X-l was discovered at X-rays by Bowyer et al. 
(1965). Early dynamical studies of the compact object sug- 
gested the presence of an accreting black hole (e.g. Bolton 
1972). The distance to Cygnus X-l is currently estimated 
to be 1.86±o J? kpc (Rcid et al. 2011). This value was de- 
termined via trigonometric parallax using the Very Long 
Baseline Array (VLB A). At this distance, the mass of the 
black hole is (Orosz et al. 2011) M = 14.8M Q . This is the 
value adopted in all calculations presented in the previous 
sections. 

The source has been observed in both a low-hard state, 
dominated by the emission of a hot corona (e.g., Dove et 
al. 1997; Gierlinski et al. 1997; Poutanen 1998), and a high- 
soft state, dominated by the accretion disk, which in this 
state goes all the way down to the last stable orbit. In the 
low-hard state, where the source spends most of the time, 
a steady, non-thermal jet is observed (Stirling et al. 2001). 
The jet is absent in the thermal state. Therefore, in this 
latter state a clearer X-ray spectrum can be obtained. 

The accretion rate and the spin parameter of the hole 
are ~ 0.472 x 10 19 g s _1 and 0.99, respectively, accord- 
ing to estimates from a Kerr plus black-body disk model 
(Gou et al. 2011). These GR models yield a spectral energy 
distribution with a maximum at E max ~ 1.6 keV. On the 
contrary, /(i?)-models with negative curvature correspond 
to a low maximum temperature, lower even than what is 
expected for the (unrealistic) case of a Schwarzschild black 
hole. Therefore, we can presume that a fit of f(R)-Kerr mod- 
els to the data would also prefer high values of maximum 
temperature, i.e., ones with non- negative curvature. Models 
with accretion rates and spin close to those obtained by 
Gou et al. (2011) and small positive curvature seem viable, 
something that is consistent with an asymptotic behaviour 
corresponding to a dc Sitter space-time endowed with a 
small and positive value of the cosmological constant. 

Deep X-ray studies with Chandra satellite might impose 
more restrictive limits, especially if independent constraints 
onto the accretion rate become available. 

7. Conclusions 

We have studied stable circular orbits and relativistic ac- 
cretion disks around Schwarzschild and Kerr black holes 
in f(R) gravity with constant Ricci curvature in the 
strong regime. We have found that stable disks can be 
formed only for curvatures in the ranges of (— oo, 10 -6 ] and 
[-1.2 x 10~ 3 ,6.67 x 10~ 4 ] in the cases of Schwarzschild 
and Kerr black holes, respectively. Current observations of 
Cygnus X-l in the soft state rule out curvature values below 
— 1.2 x 10~ 3 . Additional constrains can be imposed on spe- 



cific prescriptions of f(R) gravity. In particular, logaritmic- 
gravity prescriptions are strongly constrained by observa- 
tional data. Future high-precision determination of the pa- 
rameters of other black hole candidates can be used to im- 
pose more restrictive limits to extended theories of gravity. 
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